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Abstract 
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1 Introduction 



It was recently noticed |lj that it is possible, in the heavy mass limit, to formulate fully covariant quark models for form 
factors in which the current acts in the standard way on the heavy quark while the other quarks remain spectators. 
Following Bakamjian and Thomas (BT) ||, given the wave function in, say, the rest framed], the hadron wave functions 
arc defined in any frame through a unitary transformation, in such a way that Poincare algebra is satisfied. It was 
shown in jl| that the p 2 Isgur-Wise slope parameter was bounded in this class of model: p 2 > 0.75. 

It is known that quark models are of special value when excited states are considered, since then no other 
hadronic method is available: lattice simulation as well as QCD sum rules are practically restricted to ground states 
because they use euclidean continuation. It is therefore tempting to apply this covariant approach to the B — > D** 
decays. In f| it has been shown that these covariant quark models satisfy the Bjorken sum rule §}, §. This is not a 
trivial achievement. It comes in this class of models because the boost of the wave functions is a unitary transformation 
that keeps the closure property of the Hilbert space in all frames. 

In [| the precise formulae for the t^(w) and Ty^(w) have not been computed except for w = 1. In view of 

practical phenomenology of the B — » D** decays, Ty% i w ) an d T y\ ( w ) are obviously needed for any w. This is our 
main goal in the present paper. We do not want here to make a choice for our preferred internal wave function. These 
formulae may be used by anyone who wants to apply the BT method in the heavy mass limit to its preferred set of 
rest- frame, or infinite momentum frame or whichever frame wave functions. The Ansatz although not explicitly of 
the BT type and computed in particular frames, would probably be obtained by the BT method with the harmonic 
oscillators. The same comment applies to the calculation of B — > D** Q. In deriving the formulae for Ty\{w) and 

t^ 2 {w), we will also try to repeat, as much in a transparent way as possible, the content of the BT method applied 

to the heavy quark limit. This will be done in section |^. In section |^ we will compute the t)j\{w) and t^J 2 (w). In 
section ^ we will derive directly from the latter formulae the Bjorken-Isgur-Wise sum rule. 

2 Covariant quark models of form factors in the heavy mass limit 
2.1 Framework 

The main purpose of this model is to provide a way to implement covariance in the calculation of form factors. These 
form factors appear through matrix elements relating the initial and final states of the hadrons. We assume here a 
spectator quark model, that is, of all the quarks building the hadron, only one (labeled as 1) is the active particle. 
With this hypothesis, the amplitude of any process can be written as: 

J (2tt) (2tt) (2tt) 

where we have used the so-called one particle variables, momentum pi and spin Si, and where 0(p[,pi) s i Sl is the 
matrix element of the free one-particle operator O between one particle states of the form \p, s) (we only take into 
account two quarks since we are dealing with mesons; the generalization to n particles is straightforward 

Now, the covariance is introduced by expressing the relativistic 2-particle bound states ^ r Sl ,s 2 (Pi>P2) as a 
representation of the full Poincare group, following the Bakamjian-Thomas formalism. It turns out that this process 
is made easier if we change the variables that characterize the state VP, that is, if we introduce the total momentum 
P, (P = pi + P2), the internal momenta k%, k 2 , (M + k 2 = 0), and the internal spins S[, S 2 , leading to another wave 
function S2 (P, k 2 ) which is related to the previous one by the unitary transformation: 

fl^(ft,fi) = ^^4f§ E} d ( r ^aD(r 2 ) S2A ) ¥£ s ,(P,k 2 ) (2.1) 

s 1 ,s 2 

Let us explain the notation used in this last equation: 
* P° = \Jvt + mf 

2 Any starting frame can be chosen, for example the infinite momentum frame although some care is then needed. 
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M = 

-l 



• k 2 = Bjjipa 



• D(Ri) s . s > is the (s^ sj) matrix element for a spin 1/2 particle of the Wigner rotation Ri, which describes the 
connection spin<->intcrnal spin and which is given by: 

R l = B^ 1 B^ P] B ki (2.2) 

where B p is the Lorentz boost (Wjr, 0) — > p. 
Note that there is no apparent dependence upon ki owing to the relation which defines the internal momenta. 

Returning now to the Poincare group transformation laws, their generators are defined according to: 



p 




H = 


P° = V J 




s d 


J = 


-iP x — 




dP 



(2.3) 



boosts : K 



2 I ' dP 



PxS 
P° + M 



where 



{Si + S' 2 ) - ik 2 x - r 
afc 2 



is the spin of the meson 



These generators depend on the interaction through the mass operator M . Any given quark model corresponds to 
one choice of the mass operator M, in order, for example, to describe the mass spectrum at best. In this paper we 
do not want to enter into the discussion of the best choice for the latter operator. We simply need M to depend only 
on the internal variables and to be invariant by rotation in order to obey Poincare algebra. The mass M$ = -J (Sp 7 ) 2 
introduced above is not the real mass. It is just a tool to define properly the unitary change of variables (2.1). Only 
in the limit when the quark interaction would vanish, would Mo become equal to M. In such a limit, the generators 
in (2.2) would be the standard Poincare generators for the considered set of free particles. As such it would naturally 
obey Poincare algebra. In some sense, the trick used above in order to have the correct Poincare algebra even when 
interaction is present is to mimic, so to say, the behavior of free particles by the definition of p®, Mq. This does not 



mean that we make a weak interaction approximation. One may say that the change of variables (2.1) is inspired from 
the non-interacting case, but it may be applied whatever the strength of the interaction is. 

It is then possible to construct Q the wave function of the bound state of quarks moving with total momentum 
P and we get: 



*7*s 2 { P M = (27T) 3 S(p 1 +p 2 - P) <p Sl , S2 (k 2 ) 



(2.4) 



where ip Sl ,s 2 (k 2 ) is an eigenstate of M, S and S z , so that the relativistic wave function expressed with the momenta 
Pi and spins Si reads: 



ln° + r>° lk° k° 

9. U .M,P2) = \l m o Jf =o E WRi).^ D(R 2 ) S2iS , (2tt) 3 S(p 1 +p 2 - P) ^ iA {k 2 ) 
Notice that for P = 0, this last formula gives 

which just expresses the fact that tp is the rest-frame internal wave function. 
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Putting all these things together, we get the new following workable relation: 



(P'\o\P) = 



(p[° + p° 2 ){p1 + p° 2 ) VWM 



(2t:) 3 V M oMo JpJpJ y/pffi 

E E ^yj^)*^^" 1 )^^.^)^^!)]^,.!^^" 1 ^,.^.!,*^) (2-5) 



s', ,sL SI ,S2 



since p 2 = p 2 for the spectator quark. At this stage, this matrix element is not covariant in general, because the 
current operator O is not covariant with respect to the transformations ( |2.3| ) which contain the interaction. However, 
in the limit where the masses mi and m! x of the quark 1 tend to infinity, this last equation becomes fully covariant as 
shown in IfiJ. 



So, up to now, we have replaced the knowing of a relativistic wave function of a moving bound state by the 
knowing of a rest-frame wave function, which is an eigenstate of a mass operator (in this case the hamiltonian) which 
does not have to be relativistic. Of course, the physics lies in the complete determination of those rest-functions 
ip SltS2 (k 2 ), and in the expression of the current operator. 

Another point which must be emphasized is the unitary nature of the transformations which have been used. Unitarity 
ensures that, in any frame, all sets of wave functions are orthonormal and satisfy the closure relation. This will prove 
to be essential for the derivation of the Bjorken sum rule (see section [|). 



From now on, we will assume that the quark 1 has an infinite mass to enable the Heavy Quark Symmetries 
(flavor and spin symmetries) and the covariance of the matrix elements. Besides, we will denote the internal spins Si 
and no longer S^. 

2.2 The wave functions 

In the rest frame (P = 0), the generator of rotations writes: 

J = Si + 3 ■ ( 2 -6) 
In quark models, j equals s 2 + I and I = —ik 2 x d/dk 2 is the relative orbital angular momentum of the second 



quark, which is equal to 1 in the case of D** mesons. However, beyond quark models, the decomposition (2.6) is very 
useful in general in the heavy quark limit. Indeed, according to Heavy Quark Symmetry, the spin si of the heavy 
quark is conserved, i.e. it commutes with the Hamiltonian. It then results from conservation of the total angular 
momentum J that the angular momentum j is also conserved. It is then natural to label the lightest parity-even states 
(I = 1, P-wave in the quark models), according to the values of j, namely j — 1/2 and j — 3/2. It is known that 
all hadrons within j = 1/2, respectively j — 3/2, multiplet, are related by the heavy quark symmetry^]. Finally, when 
combined with the heavy quark spin s\, we get two distinct multiplcts: one with j — 1/2 and J p = + or 1 + , and 
the other with j = 3/2 and J p = 1 + or 2 + . The corresponding states will be denoted respectively, using the notation 
\jJ P >: 

U0+) |il + > ||1+) ||2+) 



2.2.1 Generic form of the <p's 

How are we going to write the rest-frame internal wave functions ? Recall that they are eigenstates of the mass operator 
M. In the model, this operator is assumed to be rotationally invariant, to depend only upon the internal variables 
and, of course, to conserve parity. So it commutes with J, which in turn commutes with the hamiltonian H . Now, 
the heavy quark spin symmetry implies the invariance with respect to si of the system, that is [H, si] = [M, si] = 0. 
As a consequence, [H,j] = = and the total Hilbert space TL of the <y9's can be factorized as a tensorial product 

of a spin space "Hg 1 , related to the spin si , and of a "spin-orbit" space 7ij which describes the internal spin s 2 of the 

light quark and the relative orbital angular momentum I. We must therefore build a representation of j = s 2 + I. 
But we know a priori nothing about [H, 1} or [H, s 2 ] so we cannot reproduce the same kind of argument used for the 

3 A11 these properties are valid in the model, and we shall use them in the next subsection. 
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decomposition of H into TLg 1 <g) Hj. Yet, we are concerned by the two lightest parity even multiplets (j — 1/2 and 
j = 3/2), with parity +1. They all correspond to the same value of the orbital momentum, 1 = 1. Indeed, due to 
parity conservation, I = 1 cannot be mixed with I = 0, 2, and due to the conservation of j, it cannot be mixed with 
I = 3, 5 which cannot produce j = 1/2, 3/2 when combined with s 2 = 1/2. Whence, the only way of constructing this 
representation of j is I C3> s 2 with 1 = 1. 

However, we must stress again that the radial part of the wave functions and the energies depend on j. The 
four I = 1, s 2 = 1/2 states are not degenerate in energy due to an Z • s 2 force which does not vanish in the heavy 
quark limit. Since the eigenvalue of I ■ s 2 is equal to + 1) — 1(1 + 1) — s 2 (s 2 + l))/2 = j(j + l)/2 — 1 — 3/8, we see 
that the four states will split as expected into two multiplets corresponding to j = 1/2 and j = 3/2 respectively. The 
spin and orbital parts of the wave functions combine through Clebsch-Gordan coefficients to build up the j = 1/2, 3/2 
eigenstates. As a consequence the 3 P± and the 1 P\ states will mix to produce the j eigenstates. 
Finally, we can factorize the wave function: 



* J M, S1 ,sM - (JJ P M, S1 , S2 )(P)^\\P\\M\\P\\ 2 ) 



(2.7) 



where 



(jJ P M, SuS2 ) (P) = £0' m 1/2 M -m\J M) X f~ m 1 m' 1/2 m-m'\j m) x?~ m ' Y? ' (p) (2.8) 

m m' 

Y l m (p) being the spherical harmonics (6, <f)\lm) and ^"(||p|| 2 ) a radial function which remains to be determined, and is 
normalized according to 



dp p z 
(2tt) 3 ~3 



l&(l|jfT =1 



The unusual writing of the radial part in (2.7) is a matter of convenience, meant to simplify further calculation. 
Moreover, x™ is the column matrix defined by: 



X 



+ 1/2 _ 



,-1/2 _ 



From now on, we will often identify the indices 1/2 with 1 and 
which are related to spin. Consequently, the xT equal S ms . 



-1/2 with 2 as a way of spotting the matrix elements 



Let us concentrate for awhile on the spin-orbital part of the wave function. The calculations are made easier 
if we deal with another composition of I, s\ and s 2 in order to get J, that is if, instead of writing s*i £g> (s*2 <g> Z), we use 
(s*i (g> S2) <8> l- With that decomposition, we obtain the so-called \ 2S+1 Pj) states, which can be related to the \j J p ) 
through the "6j" coefficients according to: 



\j jP ) = £r si+S2+ ' + 'V(2s+i)(2, + i) j s / s* -Jjr- / 



which reads in the present case 



^0 + ) 

fl + > 
I l§2 + ) 



V3 



I 1 Pi) 



'P 



l 3 P 2 



"Pi) 



V3 



'Pi) 



with 



G 



( 2S+1 Pj s1 , S2 m) (P) = J2( lmS M ~ m I J M)YT(p) (1/2 a 1/2 b \ S M-m) (x Q ) Sl 

m a,b 

= ^(lm5 M-m \ J M) {1/2 s x 1/2 s 2 | 5 M-m)Y{ n (p) 

i n 

since (x m ) s = 5 ms 

It may also be useful to remark that the Clebsch-Gordan coefficient (l/2sil/2s2|/S , m) can be written as 
the matrix element of a combination of Pauli matrices according to: 



1/2 si 1/2 s 2 \0m) = —(a 2 ) Sl , 



(1/2 si 1/2 s 2 | 1 m) = -j= [(eW -a) 
where we have used the standard basis: 



<T 2 



(2.9) 



e( 0) = e 2 



V2 



(fix i Cy) 



(The demonstration of these relations can be carried out by evaluating the right hand sides of (2.9) and then by 
comparing the results with the values of the Clebsch-Gordan coefficients). 

Finally by substituting the expressions of the Y[ n {p) in spherical coordinates, the \ 2S+1 Pj) spin-orbital wave 
functions write: 



( 3 ^o Sl)S2 ) (P) 



P cr 2 



( 3 ^ ai ,. a ,m) (P) 
{ 3 P2s US2 ,m) (P) 



i 

S1S2 \\p\\ 

where ej™ is the rest-frame polarization tensor, that is a symmetrical tensor with vanishing spur. 



= —i 



i 

s±s 2 \\p\\ 

o 2 



1 



We are now able to go one step further into the determination of the transition amplitude (P'\0\P). In the 
present case, the state |P) represents a pseudoscalar state 1 S$ and the state \P') one of the ( 2S+1 P/) states described 

above. In other words, the y SljS2 (fc)'s take the form: 

(p Sl ,s 2 (k 2 ) = -^=(a 2 ) SlS2 ip(\\k 2 \\ 2 ) 



<p' s > s >M) 



V2 



(X^2)s{s' 2 4>j(\\h\\ ) 



where 



^(2S + l Pj ) 



(a-k' 2 ) 



= < 



(e*M • k' 2 ) 

e (m) • (k' 2 A a) 



i 

71 



for the 3 Po state 
for the 1 Pi state 
for the 3 Pi state 
for the 3 P 2 state 
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Then, by substituting these expressions into ( |2.5D and using the relation a^DiRja^ = D(R 1 ) t , we get: 

I(pT+P 2)(p°i+P°2) VWK \[WK 



J (2ir 



(2ny V M oMo ^fpYpJ VMpI 



x ^Tr 
2 



t D(R'^) Ofoft) D(Ri) D(R^R' 2 )} ^{\\k' 2 ff v {\\k 2 f) (2.10) 



2.3 Switching to Dirac notation 

As already stated, it can be shown directly on ( |2.5| ) that the current matrix element is covariant in the heavy mass 
limit mi — > oo. We shall demonstrate covariance by changing to Dirac notation, in which the current matrix elements 
will eventually appear as manifestly covariant. This notation will also be very useful to handle later calculations. 
The idea is to insert the 2x2 matrices, that appeared in the precedent sections, into the 2x2 upper left block of a 
4x4 matrix which is then completed with zeros: (J J). This leads to the following formula: 



dp2 1 y/ U'o U a fc° fcrj fe'l ^'2 

(2tt)3 ^° p°p>° ^/k° + mi y/k° + m 2 + m[ ~^f+r7T 2 



X 16 



lTr{0(m 1+ A)(l + «(m2+^)(l + ^) (b./b, 1 ) K + #)} &(||^|| 2 )Vllfe|| 3 ) (2.11) 
where u and u' are defined by 

M u = pi + p 2 {u 2 = 1) M' u' = pi + p 2 (u' 2 = 1) 



The derivation of (2.11) is quite straightforward: it stems from 



• the expressions of the Wigner rotation matrices as a product of three boosts, see equation (2.2) 

• the expression of the matrix 0(p',p) as [Q 

0(Pi,Px) = , — —B,OB Pl 

2 2 

In this last relation, the boosts B p are written in the Dirac representation according to 

m + jfry 

p = — / 

y'2m(p + m) 

• the properties of £?„, that is, for instance: 

B u (l,0) = u (idem for B u i) 

p 2 = B u k 2 = B u > k' 2 

• the following formula, which will also be used throughout the rest of this paper in order to evaluate the terms 

B u ( lllX »)BZ l = l^{B u xf (2.12) 
(the 7^ are forming a 4- vector), and which gives for example: 

B u ~f B- x = B u (l, 0V B- 1 = f.B„(l,S) = i 

• the fact that the Wigner rotation matrices, the x's and 7 are all block-diagonal matrices, implying that they 
commute with the matrix 7 which has scalar blocks. 
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2.4 The heavy-mass limit 



Let us consider now the heavy-mass limit of (2.11). By doing this, we mean that we take m 1; 
the same time, we keep fixed the following ratios: 



m'i 



and, at 



Therefore, we have: 



EL 
~M~' 



P 

Ti 



pi 

mi 

u - 

K 

mi 
k 2 - 



P 1 
m'i 
v! - 

*X 

m'i 



As a consequence, since the scalar product of 4-vectors is invariant, we get: 

(B-W=P3.« (B- 1 p 2 )°=p 2 .v' 



We deduce from these formulae the invariance of \\B V 1 p2\\ 2 and || -B^/p2 1 



\\B- l p 2 \\ 2 = { P2 .vf - ml \\B-}p 2 f = (p 2 .v'f - ml 
Finally, the heavy mass limit of the relation ( [2.11 ) reads: 



(2.13) 



( 2S+l Pj\O\ l S ) 



1 1 



dp 2 1 



y / (p 2 .v')(p 2 .v) 



8 y/ v ° v 'o J ( 2?r ) 3 Pz \/(P2-v' + m 2 )(p 2 .v + m 2 ) 



Tr{0(l + ^)(m 2 +^) (B^B-, 1 ) (1 + f)} MWBjPlf) Vl^W) (2-14) 



( 2.14| ) is the starting point for the calculation of the transition amplitudes from a pseudoscalar state towards the 
( 2S+1 Pj) . states (£)** states): that is the purpose of the next section. 



It is important here to stress that the expressi on ( 2.14 ) is covariant. The integration measure dp 2 /p2 is 
invariant; so is the trace as obvious from expressions (3.1) for B v iX*B~, and finally the arguments of the wave 
functions are invariant from (2.13). 



3 Transition amplitudes and Isgur-Wise scaling 
3.1 Preliminary calculation 

Before calculating the spur which appears in (2.14), we must evaluate the B v i^B~} terms, for each X( 2S +i P )• ^hc 
procedure is quite automatic: 

1. We have to write each x, in the rest frame, in a covariant way, that is using 4-vectors instead of 3- vectors. That 
is realized by introducing the 4- vector = (1,0) and the expression of the Pauli matrices in term of the Dirac 
matrices: 

• ( 3 -Pq) : Starting from x (3 \ , we get the 4x4 matrix in the following way: 



L ( 3 -fo)' 

(off) 



A 
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where the 1+ 2 10 factor in the right hand side gives the form (jj) to the X- By introducing the n M mentioned 
above, this can also be written as: 



L ( 3 n>) 



1 

71 



75 



-k" 2 + fc'°7o 



l + 7o 



1 

71 



75 



-^2 + k' 2 "/fj. n>1 



l + 7o 



that is, finally, 



V3 



75 



l+7o 



• O-Pi): Introducing the 4-vector e — (0, e) (we drop the polarization index m for the time being) such that 
e • n = 0, it is not difficult to obtain: 



X 1 



• ( 3 -Pi): Another manipulation is performed here in order to deal with the cross product, but the idea is still 



the same: to add a fourth component. 

e ijk e*k' 



e ■ (k' 2 A (?) 



l3k e*k' l5lkl0 ~> -e^n^fe' 7s7(7 7 



Then, we get 



% ^ vpa n e*k' -y ~x 1 + 7 ° 



V2 



•( -Fb); I n this last case, we generalize the polarization tensor ey into which is symmetrical and has a 
vanishing spur and obeys n^e^y — 0. Thanks to this last property, the x writes: 

yt - mWp* 1+70 
^(3p 2) - 757 k 2 e p „ - 

2. We now must evaluate the expressions B V 'X^ B~, . That is achieved by inserting the factor By/B^, 1 at the right 
places in the previous formulae and extensively using ( 2.1 2| ) . Moreover, we introduce the following notations: 



B v t G-tii/ 



As a result, the By^x^B^ 1 read: 



B vx] 3 „ B j = [fc - (pa • v 1 ) ft] 75 



7 3 ^o) 



V3 
■(P2-e" 



1 + f 



--^^,V^7.75^ 



(3.1) 



1+ // 

Notice the — term in the previous relations; when contracted with the factor (1 + ft) of the spur in (2.14), 

the result will be (1 + ft) since ft 2 = 1. So, from now on, we will drop it from the expressions B v i~x^ B~} . 

3.2 The transition amplitudes 
3.2.1 Definitions 



On inserting the expressions (3.1) from the previous section into (2.14), three different kinds of integrals appear using 
the covariance of the matrix elements, ( [2.14 ). Therefore, we may introduce the following definitions regarding the 
result of these integrals: 



A U) 



dp 2 1 



y / (p 2 ■ V')(p 2 ■ V) 



1„ l|2\ 



( 27r ) 3 V° 2 y / {p 2 -v' + m 2 )(p 2 -v + m 2 ) 



B^p 2 rr<p(\\B- i P2 
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B (j) v » + B ' U) v ni 
D^n" v' v + v v v" 1 ) 



dp 2 1 



y / (p 2 ■ V')(p 2 ■ V) 



B-} P2 f)\{\\B-^r)p^ 



( 27r ) 3 P°2 y/fa ■«'+m 2 )(p 2 • v + m 2 ) 



dp 2 1 



\/(P2 ' V')(p 2 ■ V) 



(WB-^ff^WB-^p^p^ 



( 27r ) 3 P°2 y/iP2 -v' + m 2 )(p2 -v + m2) 



where the A, B's and the D's arc function of w — v ■ v' . 



3.2.2 Results 



We now have all the pieces to reduce ( J2.14 ) further for all the states (jJ p )- Straightforward calculations and a little 
diracology lead to: 



■o+ioi 1 ^) 



1 1 



-7= ^=Tr{0(l + ^) (1-^)75} 
8V3 Vvov'o 

(1 + v ■ v ) 



3D 



(1/2) 



ii + iors ) = 



1 1 



x [Tr{0(l+^)(1+^')} (ve*) + i «'„ e* v„ Tr {O (1 + 7<r (1 - 75}] 



1 1 



(1/2) 



8^ V^o W'o 4 

x [Tr{0 (1 + ^)^(1 + ^')} + i e^ CT «'^*Tr{0 (l + ^) 7 , 7ff (l-^)75}] 



l + |0|^o) 



1 1 



4^6 V^o u 'o 



= (m 2 B^+D[ 



(3/2) _,_ £,(3/2) 



-Tr{0 (1+^(1 + (we*) + 5 e^<V„ e> p Tr {O (1 + 7 . (1 - f) 75} 



4^6 \Jv Q V' 



D 



(3/2) 



-Tr { (l + f)f(l + f)} + _ e ^^^Tr{0(l+^)7p 7.(1-^)75} 



■2+\O\ 1 S ) 



■-^= -=L= (m 2 5 (3/2) + £,(3/2) + £)(3/2)\ Tr {0 (1 + ^ y, (1 _ ^ ^ , 
8V3 \Jv v' V / 



1 1 



= £>i 3/2) Tr{0 (1 + ^)7" 7" (l-^')75} 



8V3 Vvo v'o 
3.2.3 Vector and axial matrix elements 

Up to now, we did not take into account an explicit form for the current operator O. Let us see what becomes of the 
preceding relations when we use O = V\ = 7a and O = A\ = 7a 75- 



(i 0+| 14 l^o) =0 



(4 0+1 A X |%) = --i= (w A - {(1 + w t,') (m 2 + < 2 > + £# /2) ) + 3 ^ 1/2) } 

2V3 V^^o L v 7 J 



(ll+l^l^o) 



2^3 
1 1 
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_ [( W . C *K + (l-«V)eJ] 
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x{(l + vv')(m 2 B^+D[ 1/2) +Di 1/2) ) + 3£>f /2) } 

1+| A, \' So) = —L= -±= e XapT v° e*? v' T {(1 + v ■ v') (m 2 + D< 1/2) + D ( 2 1/2) ) + 3 D^X 

V12 y/v v' \ 7 J 



for the j = 1/2 multiplct, and 

(| 1+| V X l^o) = [{v ■ e*) ■ t/ - 2) - 3„a} + e*(l - t, • W ')(l + v ■ v')] 

2V6 y/v v' 



1+UaI 1 ^) 



x{m 2 B^+D^ 2) +D { 2 3/2) } 

~ -^=f (1 + » • »') eA.pr « ff ^ «' T |m 2 i?( 3 / 2 ) + D< 8 /2) + ^ 3/2) } 
2V6 V^o^o 

(| 2+| Va l^o) = j e AQ „ e*^ «- ^ |m 2 B< 3 / 2 ) + M 3 / 2 > + D^} 

<|2+|A A | 1 5o) = i^l= [(l+t/.tO^t^ - v»v»e%v\] {m 2 B^ + D { 3/2) + D^} 

2 V V o V o 1 J 

for the j = 3/2 multiplet. As expected, all these matrix elements are manifestly covariant. 
3.3 Isgur-Wise scaling 

Looking at the results of the last section, it is fairly obvious that all the matrix elements can be written using two 
differents quantities, t\i%(v ■ v') and r 3 / 2 (u • v'), following the notation introduced by Isgur and Wise in M. These 
Tj(v • u')'s are: 

r 1/2 (w') = ^= {(! + «• ,')(- 2 5 (1/2) +< 2) +^ 1/2) ) + 3^ 1/2) } 

r 3/2 (vv>) = ±{m 2 B^+D^+D^} (3-2) 

Two remarks though: our definition of the state |± 1+) differs from the one given in § by an overall minus sign; 
and, secondly, we are using another normalization of the states. As a consequence, there is an overall factor of 
mB** a/ v v' q between our transition amplitudes and the transition amplitudes written in M. 

4 The Bjorken sum rule 

Generally speaking, the Bjorken sum rule is a way of expressing duality (between quarks and hadrons), which, in the 
case of the heavy quark mass limit, is an exact duality. Mathematically, this rule reads H: 



h^(v,v')= ^J2<n,P'\O^\0,P> x{<n,P'\O v \0,P>)^J 



where mi — > oo 

= \ E [u s[ (v')0"u Sl (v)] [u^O^M* = ^ cquark (4-1) 

where hy. v (v,v') is the "hadronic tensor" which describes the transition between an initial meson state to all possible 
meson states containing a heavy quark, and where hz^ e i uark is the analog for free quarks. Note that n is used to label 
the vectors of the complete basis {|n, P >}. 

By using the expressions of the current matrix elements which are involved in the definition of h^ v {v, v') and 
expanding around w = v • v' = 1, we get the new following Bjorken sum rule: 



£(|r(;)(l)| 2 + 2|rW(l)| 2 ) = p 2 \ (4.2) 
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where p 2 is the slope of the elastic ground state Isgur-Wise scaling function and where the superscripts (n) characterize 
the radial excitations of the P wave states, to which the results of the preceding section obviously apply (these (n) 
correspond exactly to the labels n used in ( |4 .!] )). In the physical meaning and a detailed demonstration of that sum 
rule have been given. In the following, we are just going to verify that (4.2) still holds with the covariant formalism. 



4.1 The t,-(1)'s 



Checking ( |4.2Q begins with evaluating the Tj's at v ■ v' = 1. Since they are given in a covariant way, let's choose the 
particular frame of reference where: 



V 
v 







with v small 

Therefore, to the second order in v, we get: 

p ■ v' = po 
p-v ~ po 

vv' ~ 14 



(then v' = 1) 
(then v — 1 + - v 2 ) 



•r 2 



and also 



\B-hf 



= p 



p - 2p p-v 



Then, we substitute these expressions into (3.2) and expand to the first order in v the <fi functions of (3.2). Owing to 
rotational invariance, the terms of the form v/v 2 vanish and, when we take v equal to zero, we are left with: 



.(«) 



(1) 



p 2 dp ,*{n) ( 2 



(2tt)2 



(4.3) 



with p = \\p\\ and where 



F 1/2 {P 2 ) 
F 3/2 (P 2 ) 



1 

3V3 
1 

3V5 



V(P 2 ) 



V(P 2 ) 



P 



m + po 



rn 
Po 



A d( P I 2n 2 

dp 



P 



m + po Po 



A dc fi I 2^ 2 
4tt(p )PoP 
dp z 



4.2 Connection between the Fj's and the sum of the r^'s 

Let us work in the part Tij of the total Hilbert space Ti. describing the wave functions, that is, let us fix the spin of 
the heavy quark (which we know is irrelevant). Then, the partial wave function writes: 

"(p) = ^' l)m (P)( 1 m 1/2 M - m | j /x) X M " m 



In 7ij, the closure relation reads: 



^^(fltfW'V) = (2irf6(p-p f ) (4.4) 



If we multiply (4.4) by J df2 dQ! Y™ 1 (fi) Y^™ 1 (O') and sandwich the result between x" 011 the left and \ v ^ on the 
right, we obtain: 

£ y 0j (n) b 2 ) (P*) <W ( 1 m x 1/2 !/ I j M > < 1 mi 1/2 ✓ | j // ) = (2^) 3 fc^l 5 mi m , *„ „ 
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where p = \\p\\ and p' = \\p'\\. Finally, by multiplying this last relation by (1 mi 1/2 v \ J M) (1 m\ 1/2 v 1 \ J' M' ) 
and summing the result over mi, m' l7 i>, v 1 , we get: 



Ec"'(p 2 )C(p' 2 ) = 6- : 



<j(p-p') 



PF 



for each value of j. We now have all the pieces to calculate X)l r j| 2 ! starting from (4.3) and summing its squared norm, 
we get: 



<W-(p 2 )I 2 



4.3 Checking the Bjorken sum rule for p 2 



We now are able to evaluate the left hand side of (h2). After integrating by parts (the integral of the divergence 
vanishes), we get: 



5>(;)(i)| 2 + 2 2>#(i)| s 



dp 

3 j (2^y 



2n* 



tp(j3 )<p(p ) 



1 P Po + 2m o 

— —9 — p 

4 Po Po + m 



d( P ( 2\ dl P , 2\* a 2 A 
+ d^ {P) dpt (P) 4P °\ 



Regarding the right hand side of (4.2), we start from the expression of the slope p 2 given by the equation (29) of [Q: 



9 = 3 



dp 
(2^)3 



[v p V(p)]*-[VpV(^] 



dp .2 1 m 2 1 m 



(2tt) 3 l 3 4 (p ) 2 3 p°+m 



<f(p)*<P(p) 



Then, we use 



V(pa(p) = PoVip + tpVpo = 2p p^ + —if 



dp 2 



Po 



and, after another integration by parts, we get 



1 



dp 



lf(p 2 )if(p 2 



1 p po + 2m 2 

— —9 — p 

4 Po Po + m 



d f ( 2n d <P , 2x* A 2 4 



showing that (4.2) is indeed verified. 



5 Conclusion 

In this paper, we have considered the B — > D** type transitions quark models a la BT in the heavy quark limit. 
As already shown for the B — > transitions Q], these models verify covariance and heavy quark symmetry. 

Consequently, all the hadronic matrix elements between a pseudoscalar state and a P-wave meson state, both containing 
a quark with infinite mass, can be expressed using two different functions, namely t 1/2 and r 3/2 , which we have computed 
in these models: 



1 f dp 2 1 y/(p2 ■ v')(p 2 ■ v) , 



2v 3 J 1 J P2 v (P2 •« + m 2 )(p 2 • v + m 2 ) 

(P2 • v)(p2 ■ v' + m 2 ) — (p2 • «')(P2 • u' + (f • v')m2) + {l — v-v')m\ 



and 



/ \ 1 1 /" ^P2 1 y/(p2-v')(j>2-v) , ,, ,v 2 2n* // N2 2n 

= — = ; -TT7 / 7- . „ =C>3 ((jP 2 .f ) — % 0?((»2-f) — Tl 9 ) 

V; V3l-(«'«') 2 i (2tt)3 p° V(P2-«' + m 2 )(p 2 - t ; + m 2 )^ 1 ^ 7 a^VU/* , 2 , 
^— 7 [P2 • (« + v')] 2 _ (p2 ■ v)(2p 2 ■ v' ~ m 2 ) - {p2 ■ v')[p 2 ■ v' + (v ■ v')m2] - — - 
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In the two last equations the apparent poles at v ■ v' = 1 are canceled by zeros in the numerators. 

We have also checked the validity of the Bjorken sum rule already demonstrated in ||. These results are general 
for all quark models a la BT, independently of the precise dynamics. To complete the calculation, it is necessary to 
have an explicit form for the radial functions <fij and ip. Therefore, what remains to be done is to take a hamiltonian 
and then use it to solve a rest-frame Schrodinger-type equation, leading to wave functions from which the radial parts 
can be extracted. That will be done in a forthcoming paper. Of course, all the difficulties lie in the choice of the 
hamiltonian describing the meson (as an example, the resulting values of p 2 must not be too highj^]). When this is 
achieved, we will be able to compute, consistently with the Bjorken sum rule, p 2 and the rate of the semileptonic 
decay B D** Iv. 
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4 The experimental value of jp = 0.87 ± 0.12 ± 0.08 from CLEO II is not easy to compare to p 2 in our model since 1/mi corrections and 
radiative corrections should be considered. Still we would like to avoid a too large discrepancy. 
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